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Abstract 

We include nonfactorizable soft gluon corrections into the perturbative 
QCD formalism for exclusive nonleptonic heavy meson decays, which com¬ 
bines factorization theorems and effective held theory. These corrections are 
classified according to their color structures, and exponentiated separately to 
complete the Sudakov resummation up to next-to-leading logarithms. The 
nonfactorizable contributions in nonleptonic decays are clearly identified in 
our formalism, and found to be positive for bottom decays and negative 
for charm decays. Our analysis confirms that the large-A c approximaton is 
applicable to charm decays, but not to bottom decays, consistent with the 
phenomenological implications of experimental data. The comparision of our 
predictions with those from QCD sum rules is also made. 
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I. INTRODUCTION 


Nonleptonic heavy meson decays are more difficult to analyze compared 
to semileptonic decays, because they involve complicated strong interactions. 
A naive perturbative QCD (PQCD) formalism [|lj], considering dynamics be¬ 
tween the large typical scale t of decay processes, which is of order of the 
heavy meson mass, and the hadronic scale of order Aqcd, is appropriate only 
for semileptonic decays. For nonleptonic decays, a more sophisticated for¬ 
malism must be developed, which further includes dynamics between the W 
boson mass M w and the typical scale t. Recently, such a modified PQCD 
framework has been proposed 0, in which nonleptonic decay rates are fac¬ 
torized into the convolution of a “harder” function characterized by M\y, a 
hard subamplitude by t , and meson wave functions by the hadronic scale. 
The renormalization-group (RG) method and the resummation technique 
[||] are then applied to organize large single logarithms In (Mw/t) and dou¬ 
ble logarithms Iii^/Aqcd) contained in the perturbative expansions of these 
convolution functions into Wilson coefficients || and Sudakov factors OB. 
respectively. In the previous analysis of |J we have considered nonfactoriz- 
able contributions to the hard subamplitude. However, nonfactorizable soft 
gluon effects, which cannot be absorbed into meson wave functions, were ne¬ 
glected because of their complicated color flows. In this paper we shall take 
into account these soft corrections, which produce single logarithms, system¬ 
atically in the modified PQCD formalism, and improve the accuracy of the 
resummation up to next-to-leading logarithms. 

The simplest and most widely adopted approach to exclusive nonleptonic 
heavy meson decays is the Bauer-Stech-Wirbel (BSW) model |7J based on 
the factorization hypothesis, in which decay rates are expressed in terms of 
various hadronic transition form factors multiplied by some coefficients. The 
coefficients of the form factors corresponding to external W boson emission 
and to internal W boson emission are Gq = C\ + c 2 /N c and a 2 = c 2 + Ci/N c , 
respectively, N c = 3 being the number of colors. Here c t are the Wilson 
coefficients appearing in the effective Hamiltonian, 

4 Q 

H c g = —/=-V cb V* d [c 1 {n)O l + c 2 (fj)0 2 ] , (1) 

which satisfy the matching conditions Ci(Mw) = 1 and c 2 (Mw ) = 0. 0\ — 
(ci7 At &L)(^L7 /i 'Wi) and 0 2 = are the local four-fermion 


2 



operators associated with bottom decays. The form factors may be related 
to each other by heavy quark symmetry, and be modelled by different ansatz. 
The nonfactorizable contributions which cannot be expressed in terms of 
hadronic form factors, and the nonspectator contributions from W boson 
exchanges are neglected. In this way the BSW method avoids complicated 
QCD dynamics. 

Though the BSW model is simple and gives predictions in fair agreement 
with experimental data, it encounters several difficulties. It has been ob¬ 
served that the naive factorization is incompatible with experimental data 
for color-suppressed decay modes [0. Therefore, the large N c limit of < 2 ^ 2 , ie. 
the choice cq = Ci(M c ) ~ 1.26 and cq = C2(M c ) m —0.52, with M c the c quark 
mass, must be employed in order to explain the data of charm decays [j!]]. 
However, the same limit of cq = c\(M b ) ~ 1.12 and a 2 = c 2 (M b ) ~ —0.26, 
M b being the b quark mass, does not apply to the bottom case. Even after 
including the C\^/N c term such that cq = 1.03 and a 2 = 0.11, the BSW 
predictions are still insufficient to match the data. To overcome this dif¬ 
ficulty, parameters y, denoting the nonfactorizable contributions which are 
suppressed in the large N c limit 0, have been introduced 0. They lead to 
the effective coefficients 

° e f = Ci + c 2 + \i j ) a 2 S = C2 + Ci + X2^j • ( 2 ) 


The parameters x should be negative for charm decays, canceling the color- 
suppressed term 1/N C , and be positive for bottom decays in order to account 
for the observed constructive interference in B —» D^tt decays. A phe¬ 
nomenological extraction H from the CLEO data |L0| gave 


Xi (D — > K 7r) ~ x 2 (D —> Ktc) ~ —0.36 , 

Xi(B —>■ Dtt) ~ 0.05 , X2 (B —> Dtt) ~ 0.11 . (3) 


We shall demonstrate that our predictions for the nonfactorizable contribu¬ 
tions are consistent with the above values of x- 

The rule of discarding the 1/N C corrections 0, i.e. the large N c ap- 
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proximation, found its dynamical origin in the analyses based on QCD sum 
: x f° r charm decays are indeed negative, and cancel the term 1/N C 
However, the sum rule analyses also predicted negative x f° r H ie 
D + tc~ in a certain kinematic limit 0, and are thus 


roughly, 
bottom decay B° 
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in conflict with the phenomenological implications. Hence, the mechanism 
responsible for the sign change has not been understood completely, and re¬ 
mains a challenging subject in weak nonleptonic decays JT3]. In our formal¬ 
ism with the soft gluons taken into account, the nonfactorizablc contributions 
can be clearly identified. We shall investigate how the soft corrections modify 
the previous predictions for heavy meson decays @. ||, and try to explore the 
dynamical origin for the sign change of the nonfactorizable contributions in 
bottom and charm decays in a unified viewpoint. 

In Sect. II we briefly review the derivation of the three-scale PQCD fac¬ 
torization formulas for heavy meson decays, and include the nonfactorizable 
soft gluon effects into the formulas. The numerical results, along with a de¬ 
tailed comparision with those from the phenomenological and QCD sum rule 
analyses, are presented in Sect. III. Section IV is our conclusion. 


II. SOFT GLUON CORRECTIONS 


The construction of the modified PQCD formalism is as follows. Before 
including QCD corrections, nonleptonic heavy meson decays are described 
by a full Hamiltonian of four-quark current-current operators, such as 


4 Gf 

H = ^-V cb V* d {c L ^^b L ){d L Yu L ) 


(4) 


for the B —*■ D W 7 r(p) decays. Consider one-loop corrections to a tree-level 
heavy quark decay amplitude without integrating out the W boson. The am¬ 
plitude is ultraviolet finite because of the current conservation (a conserved 
current is not renormalized) and the presence of the W boson line. However, 
these corrections give rise to infrared divergences, when the radiative gluons 
are soft or collinear to the involved light quarks. We classify the higher-order 
diagrams into the reducible and irreducible types ||I4[: The former contains 


double logarithms from the overlap of soft and collinear divergences, while 
the latter contains only single soft logarithms. 

We factorize the infrared sensitive contributions according to Fig. 1(a) 
first, which are collected by the appropriate eikonal approximation for quark 
propagators. The diagram in the second parentheses is absorbed into a meson 
wave function 0(6, /i), if it is two-particle reducible, or into a soft function 
1/(6, n), if it is two-particle irreducible as exemplified by Fig. 1(a). Both 0 
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and U depend on a renormalization scale //, since the eikonal approximation 
brings in ultraviolet divergences. The variable 6 , conjugate to the transverse 
momentum kr carried by a valence quark, can be regarded as the spatial 
extent of the meson. 1/6 is then the hadronic scale, which will play the role 
of an infrared cutoff for loop integrals below. 

The diagrams in the first parentheses of Fig. 1(a) involve scales above 
1/6, and need further factorization. We express the full 0(a s ) diagram into 
two terms, with the first term obtained by shrinking the W boson line into a 
point, and the second term being the difference between the full diagram and 
the first term. Evidently, the former is characterized by the scale t -C M\y 
introduced before, and the latter by momenta of order My y. We factorize the 
contributions with the scale M\y according to Fig. 1(b), grouping it into a 
“harder” function H r (Mw, /x) (not an amplitude). The remaining diagrams 
in the last parentheses, because of operator mixing, correspond to the local 
four-fermion operators 0\ or 0 2 in Eq. ( 0 )- They are absorbed into a hard 
decay subamplitude H(t, /i). Similarly, shrinking the W boson line brings in 
ultraviolet divergences, and thus both H r and H acquire a /i dependence. 

Therefore, we arrive at the factorization formula for the decay amplitude 


M = H r (M w ,n) 0 H(t,n ) 0 0(6, fi) 0 U(b,fjb ) , 


(5) 


where 0 represents a convolution relation, since t and 6 will be integrated 
out at last. Assume that 7 n r , "}</>, and 7 jj are the anomalous dimensions of 
H r , 0, and U, respectively. The anomalous dimension of H is then 7 h = 
— {lH r + 7 c/> + 7 u)i because a decay amplitude is ultraviolet finite as stated 
above. A RG treatment of Eq. (|) leads to 


M = H r (Mw, M w ) 0 H(t, t) 0 0(6, 1 / 6 ) 0 U(b, 1/6) 




which becomes explicitly fi- independent. The two exponentials, describing 
the two-stage evolutions from Mw to t and from t to 1 / 6 , are the consequence 
of the summation of large In (Mw/t) and ln(f 6 ). The first exponential can 
be easily identified as the Wilson coefficient c(t). Note that the Wilson 
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coefficient appears as a convolution factor here, instead of a multiplicative 
factor in the effective Hamiltonian, Eq. (0). In the conventional approach 
to exclusive nonleptonic heavy meson decays [jTJ, which is based on Eq. 
a value of /j must be assigned, and thus an ambiguity is introduced [|1^ 


A merit of our formalism is then that it does not suffer the scale setting 
ambiguity. 

Equation (|6j) sums only the single logarithms. In fact, there exist the 
double logarithms ln 2 (P + 6) in the meson wave function 0, which arise from 
the overlap of collinear and soft divergences, P + being the largest light-cone 
component of the meson momentum. Hence, 0(6, /i) in Eq. (j5|) should be 
replaced by 

0(P + , 6, n) = 0(6, n) exp[—s(P + , 6)] , (7) 

where the exponential e~ s , the so-called Sudakov form factor, comes from the 
resummation of the double logarithms. For the detailed derivation of Eq. © 
and the complete formula of s, refer to 0 §• Below we shall approximate 
H r (M Wl M w ) by its lowest-order expression H= 1, and evaluate H{t,t ) 
perturbatively, since all the large logarithms have been grouped into the ex¬ 
ponents. For simplicity, we set the nonperturbative initial condition U (6,1/6) 
of the RG evolution to unity due to the strong suppression at large 6 |Hj, 
and neglect the 6 dependence in another nonperturbative initial condition 
0 ( 6 , 1 / 6 ). 

At the tree level, six types of diagrams contribute to the hard decay 
subamplitude of B~{P\) —> D° (P 2 )7t“ (P 3 ) as shown in Fig. 2, where the 
momenta are chosen as 


Mi 


P 1 = ^|(1,1,0), P 2 = ^|(l,r 2 ,0), P 3 = ^|(0,l-r 2 ,0), (8) 

with r = Mr,/M b , Mb {Mb) being the B (. D ) meson mass. These dia¬ 
grams represent external W emissions, if the four quark operators are 0\ in 
Figs. 2(a) and 2(b), and 0 2 in Figs. 2(c)-2(f). They are internal W emis¬ 
sions, if the four quark opeartors are 0 2 in Figs. 2(a) and 2(b), and 0\ in 
Figs. 2(c)-2(f). At this level, only Figs. 2(e) and 2(f) give nonfactorizable 
contributions. We shall argue that after including the soft gluons, Figs. 2(c) 
and 2(d) become nonfactorizable. 

We discuss the color structure of soft gluon exchanges. If a soft gluon 
crosses the hard gluon vertex, the color structure is given by 

(^“Wi ( T %R 2 ( T< 1 )a 2 a ' 2 (T H )a' 2ai 


M 


M; 


B , 
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1 


= o(r 


H\ 


)a2b2 


(T )t. 


a i 


2iV r 




a2ai ? 


(9) 


where T a (T^) is the color matrix associated with the soft (hard) gluon. Con¬ 
tracted with S ai b 1 S a2 b 2 from the color-singlet initial- and final-state mesons, 
the first term, denoting an octet contribution, diminishes. The second term 
leads to a factor —l/(2iV c ) without changing the original tree-level color flow. 
If a soft gluon does not cross the hard gluon vertex, it introduces a factor 
T a T a = Cp = 4/3. It will be shown that each type of hard diagrams acquires 
different soft corrections, which must be organized separately. 

For the external W emissions, the soft function receives two contributions 
as shown in Fig. 3(a): 

u [e) = ~\I(puJh) + I(h,k^)}. (10) 

Employing the eikonal approximation for the quark propagators, the loop 
integral / is written as 


I(u,v) = —ig 2 [i 


d 4 e l u a v 0 N a(3 
(27r) 4_e u ■ Iv ■ l l 2 


( 11 ) 


with the tensor 


-^a/3 9a(3 


n a l /3 + npl c 
n ■ l 


+ n 


Lb 

(n • l)' z 


( 12 ) 


from the gluon propagator in the axial gauge n-A = 0, n = (1, —1, 0 ) being a 
gauge vector. The other soft contributions either vanish because of the color 
flow in Figs. 2(a) and 2(b), or cancel by pairs in Figs. 2(c)-2(f). 

For convenience, the dimensionless vectors associated with the corre¬ 
sponding meson momenta are assigned as 


Pi = (1,1,0), p 2 = (l,r 2 ,0) , 
h = (0,1,0), fc 2 = (1,0,0), (13) 


for the b quark, the c quark, the light valence quark in the B meson, and the 
light valence quark in the D meson, respectively. The dimensionless vector 
associated with the pion is then p 3 = (0,1,0). All the above kinematic 
variables can be copied to charm decays directly, such as > Kir, by 
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choosing r = Mk/Md, Mk being the kaon mass. Varying r, we can study 
how the soft corrections modify the predictions for bottom and charm decays. 
Concentrating only on the pole terms, we obtain 


CM. 1 

I(Pi,P2) = - S --[f(Pi,P 2 )-f(pi,n)-f(p 2 ,n) + l], 

7T e 


I(ki,k 2 ) = 


a. q 1 


7 r e 


, k\ ■ k 2 1 (h- n 

In-In --— 




with 


f(u,v ) = u-vl dx[x 2 u 2 + 2x(l — x)u ■ v + (1 — x) 2 v 2 ] 1 
Jo 


i 2 1 (k 2 ■ n) 2 

-n + 1 

2 n 2 


A straightforward calculation gives 

a, 1 / 4r 2 


I(Pl,P2) = 
I(ki,k 2 ) = - 


7 T 6 \ r 

Ots 1 


4 -1 


In r + 1 


S 

7r e 


from which we extract the anomalous dimension 

2 r 2 


(e) 1 

Tu = 


7r N c \ r 4 — 1 


In r + 1 


(14) 

.( 16 ) 

(16) 

(17) 

(18) 

(19) 


For the internal W emissions, the analysis is more complicated. Soft 
corrections to Fig. 2(a) and 2(b) are similar to Eq. ([TO]) but with p 2 and k 2 
replaced by p 3 , that is, 

c (il) = ~y(Pi,Ps) + /(*!,Pa)] ■ (20) 

ffere k\ should be chosen as k± = (1,0,0) to render the second term mean¬ 
ingful. Other soft corrections, again, vanish because of the color flow. To 
Figs. 2(c) and 2(e), six soft gluon exchange diagrams contribute as shown in 
Fig. 3(b). The soft function is given by 

U {i2) = -^[ I(p 1 ,p 2 )-I(pi,k 2 )+I(p 1 ,p 3 ) + I(k 1 ,p 3 )] 

+C F [I(k 2 ,p 3 ) - I(p 2 ,p 3 )] ■ (21) 
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Notice the minus signs in front of I(p \, k 2 ) and /(p 2 ,p.a)- A minus sign 
appears, when the soft gluon attaches a quark and an antiquark |T6|. If the 
D meson is massless, pair cancellations occur between /(pi,P 2 ) and /(pi, k 2 ) 
in the first brackets, and between the two terms in the second brackets. £/b 2 ) 
then reduces to U^ ll \ The soft function associated with Figs. 2(d) and 2(f) 
is 


U {a) = C P [/(p 1 ,p 2 )-/(p 1 ,fe)] 

-^-(-f(Pl,P3) + I(h,n) + I( k 2,P3 ) - -T(P2,P3)] , (22) 

according to Fig. 3(b). We present only the explicit expression of /(pi,p 3 ), 


I(Pl,P3 ) = 


1 


7 r e 


1, (Pi • Pa) 2 1 , (P3 • w) 


In „ 

2 p? 


In ■ 


rr 


/(pi,n) + 1 


.(23) 


All other /’s can be derived simply from Eqs. (14), (|15D and (^3|) with ap¬ 
propriate replacement of the kinematic variables. Following the similar pro¬ 
cedures, we have 


(il) 

Tu = 

Q( s 

i 

7T A c ’ 

(*2) 


' i 

7f/ = 

7T 


(i3) 

Q( s 

' l 

7cr = 

7r 

Ac 


2 r 2 \ 2 r 2 

lnr + 1 — Cp —-- lnr 


4 — 4 

„2 


r 2 + 1 


Ah \ r 2 + 1 


r" \ 4r 2 

In r + 1 — Cp —-r In r 


r 4 — 1 


(24) 


III. FACTORIZATION FORMULAS 

The decay rate of B~ —■> Z) 0 7r _ is [^J 
1 


r = l28 - G F\ V *\ 2 \ Vud \ 2M B 


(1 - r 2 ) 3 


\M\*, 


(25) 


with the decay amplitude 

A4 = /tt[( 1 + ?")£+ — (1 — r )£-] + + -A4 e + A4i , (26) 
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fo and f n being the D meson and pion decay constant, respectively. The 
form factors £± associated with the external W emissions from Figs. 2(a)-2(d) 
are given by 


/* 1 /* oo 

£+ = 16 nC F y/rMg / dxidx 2 / bidbib 2 db 2 cj) B (x 1 )(j) D (x 2 )a s (t e )ai(t e ) 

Jo Jo 

x[(l + C + x 2 r)h e (x 1 ,x 2 ,b 1 ,b 2 ,m e ) + (r + ( + x 1 )h e (x 2 ,xi,b 2 ,b 1 ,m e )] 
x exp [~S B (t e ) - S D (t e ) - Su\t e )} , (27) 

/• 1 /•oo 

= IQttCfVtMb / dxidx 2 / bidbib 2 db 2 (j) B (xi)(j) B (x 2 )a s (t e )ai(t e ) 

Jo Jo 

x(_[x 2 rh e (x!,x 2 , b u b 2 , m e ) - x 1 h e (x 2 , x u b 2 , b u m e )] 
x exp[-S B (te) - S D (t e ) - sffl(t e )] , 

with the constants |5|] 


(28) 


2 


3 

2 + 


IT) - 1 

r/ + 1 


^"2 


C- = 


^“2 + 




(29) 


rj = (1 + r 2 )/( 2 r) is the maximal velocity transfer involved in the decay 
process. The form factors £; = £ ?;] + £ i23 associated with the internal W 
emissions are 

r 1 roc 

€n = / dx ± dx 3 / b^b^db^Bix^^x^a^t^c^U) 

Jo Jo 


x 


[1 + £ 3(1 - r 2 )}hi{x 1 , x 3 , bi, b 3 , m*) + xir 2 hi(x 3 , xi, b 3 , 61 , 


x exp[-5' s (i i ) - ^(fj) - S^ 1 ^)] , 


(30) 


from Figs. 2(a) and 2(b), and 

/* 1 /•oo 

623 = 16 irCpy/rMl / dx 3 dx 3 / b^b^db^six^^ix^oisiti) 

Jo Jo 

x |[1 + x 3 (l - r 2 )]hi(x\, x 3 , b\, b 3: rrii) exp[-S^ 2 ) (f;)] 
+x 1 r 2 h i (x 3 ,xi, b 3 , 6i,m i )exp[-S'^ 3) (t))]} 
x exp [~S B (ti) - S V (U)] , 


Cl(tj) 

N 


( 31 ) 
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from Figs. 2(c) and 2(d). Strickly speaking, should not be classified as a 
factorizable contribution because of the irreducible soft gluons that attach 
the B and D mesons and the pion and D meson. Only at the lowest order 
of soft corrections, i.e. Su = 0, is it a factorizable contribution. We shall 
consider this point, when identifying the nonfactorizable contributions below. 
The exponents 

Sb(/j) = s(x 1 P^,bi)+ 2 [ —7 (a s (/2)) , 

Jl/h /U 

Sd{p) = s(x 2 P 2 + , 62 ) + s((l ~x 2 )P 2l b 2 ) +2 f j(a 8 (fi)) , 

■h/b 2 n 

SM = s(x 3 P^,b 3 ) + s((l -x 3 )P^,b 3 ) + 2 f j(a s (fi)) , (32) 

Jl/b 3 fl 

correspond to the summation of the reducible radiative corrections grouped 
into the B meson wave function </>#, the D meson wave function (f>£>, and 
the pion wave function (f) n , respectively. The quark anomalous dimension 
7 = — oc s /tt, is related to 7 ^ = 27 introduced before. The exponents 

Su\v) = f , (33) 

J W fl 

for j = e,il,i2, and i3, correspond to the summation of the irreducible 
corrections grouped into the soft function U, where the lower bound w = 
min(l/&j) is chosen to collect the largest soft logarithms. The wave functions 
satisfy the normalization 


J o (j>B,D,n(x)dx = , (34) 

f b being the B meson decay constant. 

In Eqs. ©, o, ® and (|3T| ) the functions fi’s, the Fourier transform 
of the lowest-order H from Figs. 2(a)-2(d), are given by 


h e (x u x 2 ,b lj b 2 ,m e ) = K 0 {^/x[x^m~ e b\) 

x [9(bi - b 2 )K 0 (y/x 2 m e bi) I 0 {y/x 2 m e b 2 ) 

+6{b 2 - bi)K 0 ( y/x 2 m e b 2 ) /„ (y/x 2 m e b 1 )\, (35) 
hi{xuXzMM^mi) = h e (xi,x 3 ,bi,b 3 ,mi) , (36) 
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with m e = Mg and m* = (1 — r 2 )Mg. The hard scales i’s take the maximum 
of all energies involved in H : 


t e = max( v /xim e , ^/x 2 m e , l/b u l/b 2 ) 

ti = max( 1 /xi rrii, ^x 3 mi, l/&i, l/& 3 ) . (37) 


It is an important feature that the Sudakov form factor e -s exhibits a strong 
suppression in the large b region. Hence, Sudakov suppression guarantees 
that the main contributions arise from the large t region, where the running 
coupling constant a s (t) is small, and perturbation theory is relatively reliable. 

The factorization formulas for the nonfactorizable external and internal 
I'T-emission amplitudes M e and AT;, respectively, contain the kinematic vari¬ 
ables of all the three mesons. The integration over 63 can be performed 
trivially, leading to 63 = b\ or 63 = b 2 . Their expressions are 


M e = 32irV2NC F v / rMz f [dx\ [ bidbib 2 db 2 (p B (xi)(j) D (x 2 )(j) n (x 3 ) 

Jo Jo 

^ expf-g^) 163=6, - 4 6) (4 1} )] 

x[(l -r 2 )(l - ,r 3 ) - Xi + (r - r 2 )(x 1 - x 2 )]h£\xi, bi) 

-a s (t {2) ) C2 ^ ^ exp[-^(4 2) )|b 3= b 2 - 4 e) (4)] 

x[(2 - r)x 1 - (1 - r)x 2 - (1 - r 2 )x 3 ]/i( 2) (x;, 6*) , (38) 

Adj = 327 tv / 2 ~NC F \frM 2 B f [dx\ [ b 1 db 1 b 2 db 2 cj) B (xi)(j) D (x 2 )(j) n (x 3 ) 

Jo Jo 

x {« s (4 1) )^-^exp[-5(T- 1) )|6 3 =6 1 - 4 2) (T (1) )] 

x[xi -x 2 - x 3 (l - r 2 )]hf\xi,bi) 

+a a (4 2) ) Cl( ^ ^ exp[-5(4 2) )|6 3 =& 1 - 4 3) (4 2) )] 

x[(xi + X 2 )(1 + r 2 ) - l]hf } (x;,T)| , (39) 


with [dx\ = dx\dx 2 dx 3 and S' = S(b + ST + SV. The functions hP\ j = 1 and 
2, appearing in Eqs. fl38[) and (|39|) , are derived from Figs. 2(e) and 2(f): 

hf = [e(h - b 2 )K 0 (BMsbt) Jo (BM B b 2 ) 


12 










X 


(40) 


+e(b 2 - h)K 0 ( BM B b 2 ) Jo (BMeh)] 

( K 0 (BjM B b 2 ) for Bj> 0 \ 
l for Bj <0 ) ’ 

h\ 3) = [e{b 1 -b 2 )K 0 (DM B b 1 )I 0 {DM B b 2 ) 

+6(b 2 - bi)K 0 ( DM B b 2 ) J 0 {DMM] 

( K 0 (D 3 M R b 2 ) for Dj> 0 \ 

X V f H^dDjlMsh) for Dj <0 ) ’ 

with the variables 

B 2 = xia; 2 (l - r 2 ) , 

Bl = (xi - x 2 )x^{l - r 2 ) + xix 2 [l + r 2 ) , 

B 2 = XiX 2 (l + r 2 ) - (x x - x 2 )(l - x 3 )(l - r 2 ) : 

D 2 = XiX 3 (l — r 2 ) , 

L > 2 = (xi - x 2 )z 3 (l - r 2 ) , 

B)\ = (xi + x 2 )r 2 - (l - xi - x 2 )x 3 (l - r 2 ) . (42) 

Similarly, the scales d- 5 ') are chosen as 

tP = max(5M B , \Bj\M B , l/h, l/b 2 ) , 

tp = m.ax{DM B ,\Dj\M B ,l/bi,l/b 2 ) . (43) 


IV. RESULTS AND DISCUSSIONS 


In the evaluation of the various form factors and amplitudes, we adopt 
Gf = 1.16639 x 10 -5 GeV -2 , the decay constants f B = 200 MeV, fn = 220 


MeV, and f n = 132 MeV ||T0|j , the CKM matrix elements \V c b\ = 0.040 and 
\V u d\ = 0.974, the masses M B = 5.28 GeV and M D = 1.87 GeV ( 
meson lifetime t b - = 1.68 ps [0, the B and D meson wave functions 


the B 

m 


<I>b,d( x ) = 


Nb,d 


x(l — x)‘ 


16n 2 M 2 dM + C B:D ( 1 - x) ’ 


(44) 
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and the Chernyak-Zhitnitsky pion wave function ||20|j, 


, x 5v/6 
Ax) = ~—f n x{l 


x)(l — 2x)“ 


( 45 ) 


The normalization constants Nb = 590.8 GeV 3 and No = 92.85 GeV 3 , 
and the shape parameters Cb = —27.6 GeV 2 and Co = —3.372 GeV 2 are 
determined by fitting the predictions from Eq. (p5|) to the data of the decay 
rates of B —> E3- 

Our formalism can be applied to the decay D + —>■ K°n + directly, which 
occurs through a similar effective Hamiltonian but with the four-fermion op¬ 
erators O x = {s L 'y /1 c L )(d L 'y ti u L ) and 0 2 = (d L 7 At c L )(s L 7 /i M L ). The expression 
of the decay rate T is also similar to Eq. (^5j), but with the CKM matrix 
element \V CS \ = 1.0 substituted for V r j,\ , and the meson masses M D and 
Mk = 0.497 GeV [^l[ for Mb and Mb, respectively. In all the form factors 
and amplitudes the kinematic variables of the B ( D) meson are replaced by 
those of the D ( K ) meson. The D + meson lifetime is r^+ = 1.05 ps f2lj, and 
the D meson wave function has been defined above. For the kaon, we adopt 
the decay constant fx = 160 MeV, and the wave function |20|, 


4>k{x) = [3.0(1 — 2x ) 2 + 0.4] . 


(46) 


Because of the smaller D meson mass, the transverse degrees of freedom 
are more important in the definitions of the hard scales t. Hence, we choose 
the maximum of the scales 1/bi for the arguments t of the Wilson coefficients. 
In this case Sudakov suppression is weaker, and insufficient to diminish the 
contributions from the region with t close to Aqcd, where Ci^t) diverge. To 
have meaningful predictions, we choose t = max(l/ 6 j,t c ) in the numerical 
analysis, t c = (1 + 6 )Aqcd, such that the Wilson coefficients are frozen at 
ci, 2 (G) as max(l/ 6 j) < t c . The parameter e = 0.0000227 is determined from 
the data of the decays D —> 
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After including the soft gluons, the factorizable external IV-emission con- 
tributon an( ] the factorizable internal VE-emission contributon to the 
decay amplitude M. in Eq. (E 6 l) should be identified as 


( 47 ) 


xi f) = UK 1 + r )£+ - (t - r )£-}, 

Xi ^ = fD(K)£i\su=0 , 
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where ^( 5^=0 denotes the lowest order internal IF-emission contributions 
obtained from Figs. 2(a)-2(d) by turning off the soft corrections. The non- 
factorizable external IF-emission contributon and the nonfactorizable 

internal VF-emission contributon are then 

X ( " n = M,, 

XT''’ — fD(K)(£i ~ ■■ I+ Mi . (48) 


We shall show that our predictions of xi ,2 are positive and negative for the 
B and D meson decays, respectively, consistent with the phenomenological 
arguments in Eq. O 0-1 


. However, the positive nonfactorizable contri¬ 
bution Xi f° r bottom decays is contrary to the QCD sum rule predictions 
, [23|, [24|, which we shall comment on later. 

For r = M d !M b = 0.354, we have the anomalous dimensions Xu — 


—0.421a s /7r, 7 fj 1> = — 0 . 333 a s / 7 r, 7 ^ = —0.730a s /7r, and 7 ^ 3) = 0.413a s /7r. 
For r = Mk/Md = 0.265, we have xlf — — 0.396ct s /7r, 7 ^ = — 0 . 333 ck s / 7 t, 
7 ^ 2 ' ) = —0.628a s /7r, and 7 ^ 3) = 0.196a s /7r. The results of the various form 
factors and amplitudes for the decays B~ —> D°ti~ and D + —> K 0 tc + are 
exhibited in Table I. The rows entitled by Su 7 ^ 0, whose values match the 
data, are derived with the soft corrections taken into account. Those by 
Su = 0 only help to extract the nonfactorizable contributions X^ n ^\ an d to 
investigate the importance of the soft corrections. Because the D meson is 
lighter, the scales t can run to a lower value as indicated by Eqs. (|37|) and 
(f43|). The exponentials e S,J , which basically act as enhancing factors, then 
amplify the contributions from this region with smaller t, where the Wilson 
coefficients are larger. Therefore, the soft gluon effects are more important 
in charm decays as shown in Table I. 

The factorizable external VF-emission contributions are positive in both 
the B and D meson decays, and their magnitudes increase, after including 
the soft corrections with 7 ^ < 0. changes sign, since the Wilson coefficient 
C2 hi becomes so negative, when evolving from the characteristic scale of 
the B meson decay to that of the D meson decay, that it overcomes the pos¬ 
itive Ci/N c in £ i 2 3- Their magnitudes also increase because of 7^, 7 u^ < 0. 
The real parts of the nonfactorizable amplitudes A4 e (M-i) are always nega¬ 
tive (positive) due to the negative c 2 (positive C\/N c ) in the B and D meson 
decays. At first sight, these observations differ from the phenomenological 
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extractions in Eq. 


As argued before, the nonfactorizable contribu¬ 


tions should be appropriately identified according to Eq. (f48f), whose results 
are listed in Table I. We have Re(x | n ^) = +0.0189 GeV for the B meson 
decay and Re(xij — —0.4715 GeV for the D meson decay. It is easy to 
attribute the sign change of xF^ the stronger enhancement of £* by the 
soft corrections in charm decays. 

To compare our predictions with Eq. (H), we present in Table I the BSW 
parameters ay, a 2 , C2X1, and C1X2, corresponding to xf£\ Xi\ and Xi\ 
respectively. The signs, except for those of Re(x and C2X1 associated 
with the D meson decay, are consistent. The discrepancy between the signs 
of Re(xF^) and c 2 Xi for the D meson decay will be resolved, if two-loop soft 
functions U are considered. At this order, two soft gluons attach the valence 
quark lines of, for example, the B meson and the pion in Figs. 2(a) and 2(b), 
for which the associated color traces do not vanish. Then the form factors 
£± should be classified as being nonfactorizable, and xF^ be redefined by a 
similar expression to Eq. (f48|) . It is expected based on Table I that Re(xF^) 
for the B meson decay remains negative, while that for the D meson decay 
changes sign. 

At last, we come to the comparision of our predictions with those derived 
from QCD sum rules |T3|. The conclusion that the nonfactorizable contribu¬ 
tion Xi hi bottom decays is negative holds only in the kinematic limit Mb, 
Mb — > 00 but with AM = Mb — Mb fixed ||13|| . I 11 this limit both the B 
and D mesons are at rest, and soft corrections should be more important. 
The reason is similar to that for the larger soft effects in charm decays given 
before. For extremely heavy mesons, the wave functions peak at x —> 0 [|5|], 
and the characteristic scales t in Eq. (|37|) becomes smaller. Then the form 
factors £ may increase by a significant factor, since the contributions from 
the region with large Wilson coefficients are enhanced by the exponentials 

p-Su 


Following Eq. 


it is possible that the sign of %(”/) f or bottom 


decays is reversed. However, we argue that the above kinematic condition is 
inappropriate for realistic B meson decays, and thus the conclusion in [|T3|] is 
in doubt. 

In this paper we have not only completed the Sudakov resummation up 
to next-to-leading logarithms, but also investigated the soft gluon effects in 
heavy meson decays based on the three-scale PQCD factorization theorem. 
By identifying the nonfactorizable contributions carefully, we are able to 
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explain the variation and especially the sign changes of the amplitudes in 
bottom and charm decays. Certainly, these subjects still need more thorough 
studies in order to have a full understanding of their dynamical origin. 

We thank H.Y. Cheng for a helpful discussion, and T.W. Yeh for his as¬ 
sistance in the numerical analysis. This work was siipported by the National 
Science Council of Republic of China under the Grant No. NSC-86-2112-M- 
194-007 and NSC-86-2112-M-002-010-Y. 
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TABLE I. The various form factors and amplitudes for the decays B 
D° 7r” and D + —> K°n + . The unit is 10 -3 GeV. 


B —> Dir 

external 1+ 
(factorizable) 

fuii 

M e 

Mi 

Su = 0 

106.5 

2.5 

-5.8 + 19.8* 

18.5-12.5 * 

Su^ 0 

108.5 

2.6 

-5.8 + 20.0/ 

18.8-11.0* 

D -> ATtt 

external IE 
(factorizable) 

Mi 

Me 

Mi 

Sc/ = o 

267.0 

-21.3 

-18.8 + 19.0* 

37.8 - 36.7* 

Sc/7^0 

1075.0 

-529.3 

-18.5 + 13.6/ 

36.5 - 30.3/ 



A.i 

-^G/) 

Y (nf) 

A.i 


(«i) 

(a 2 ) 

(c 2 Xi) 

(ci+ 2 ) 

—» Dtt 

108.5 

2.5 

-5.8 + 20.0/ 

18.9-11.0* 


(1.03) 

(0.11) 

(-0.01) 

(0.12) 

D -> ATtt 

1075.0 

-21.3 

-18.5 + 13.6* 

-471.5 - 30.3/ 


(1.09) 

(-0.1) 

(0.19) 

(-0.45) 
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Figure Captions 


Fig. 1. (a) 0(a s ) factorization of infrared and hard contributions, (b) 
0(a s ) factorization into a “harder” function, a soft function and a hard 
decay subamplitude. 

Fig. 2. Lowest-order diagrams for the hard subamplitude. 

Fig. 3. Soft gluon corrections to Fig. 2. 
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Fig.3 



